Effective contact model for transport through weakly-doped graphene 
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Recent investigations address transport through ballistic charge-neutral graphene strips coupled 
to doped graphitic leads. This paper shows that identical transport properties arise when the leads 
are replaced by quantum wires. This duality between graphitic and metallic leads originates in 
the selection of modes with transverse momentum close to the K points, and can be extended to a 
wide class of contact models. Among this class we identify a simple, effective contact model, which 
provides an efficient tool to study the transport through extended weakly-doped graphitic systems. 
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C3 



-a 
o 

> 
On 
O 
(N 



O 



i 

C 

o 
o 



X 



I. INTRODUCTION 

Since the recent break-through in its fabrication; 1 -! 2 ^^ 
graphene, the atomically thin two-dimensional hexago- 
nal arrangement of carbon atoms, has caught the ex- 
citement of experimentalists and theoreticians alike be- 
cause it possesses unique electronic properties which orig- 
inate from the two conical points of its Fermi surface, 
located at the K points of the hexagonal Brillouin zone. 
In a simple model, a clean graphene sheet is described 
by a tight-binding Hamiltonian H = —"fJ2(ij) c l c ji 
in which next neighbors (ij) on the hexagonal lattice 
(with lattice constant a) are connected by a hopping 
matrix element 7. The lattice supports Bloch waves 
with wave numbers k x , k y and dispersion relation E — 

±7|1 + 2e l3fe * a/2 cos(\/3fc a a/2)|. At the Fermi energy 
E = of undoped, charge neutral graphene, the Fermi 
surface shrinks to two inequivalent points (k x ,k y ) = 
(|^-,±!^=^). In the vicinity of these points the disper- 
sion relation can be linearized, and then becomes conical 
with slope \dE/dk\ = 3ja/2. The ensuing low-energy 
theory is described by a Dirac Hamiltonian^ 

By varying the Fermi energy via a gate voltage across 
the charge-neutrality point, graphene offers the unique 
possibility to switch the doping of the system from n 
to p, where the doped charge density can be changed 
continuously. The first graphene-based field-effect tran- 
sistor was demonstrated in Ref. [l|. The mobility of care- 
fully fabricated graphene flakes already exceeds 50,000 
cm 2 /Vs. 4 Motivated by the ensuing prospect of graphitic 
electronic devices with properties very much different 
from semiconductor-based technology, a number of recent 
works have explored the phase-coherent transport prop- 
erties of finite segments of graphene, connected via leads 
to electronic reservoirs. This resulted in the discovery 
of a modified quantum Hall effect 2 ^ with a half-integer 
sequence of Hall conductance plateaus. It was further- 
more predicted that weak-localization corrections to the 
conductance depend sensitively on symmetries preserved 
or broken by the disorder £i£ Clean graphitic samples ex- 
hibit a finite conductivity of the order of the conductance 
quantumi^^iiSi and a shot noise identical to a disor- 
dered quantum wire£ Transport across p-n junctions was 



theoretically studied in Ref. Ill A similar but modified 
set of unique transport phenomena is found for graphene 
bilayers, which furthermore offer the possibility to open 
up a controllable gap by an electric field effect (see Ref. 
Il2l and references therein) . 

A problem untouched by the recent transport investi- 
gations is the role of the leads connecting the graphitic 
sample to the electronic reservoirs. Some of the recent 
theoretical transport studies model the leads as strips 
of doped graphene. This type of lead supports distinc- 
tively different sets of modes than a conventional quan- 
tum wire — for the same transverse mode profile, propa- 
gating modes in a conventional quantum wire often cor- 
respond to evanescent modes in a graphene strip, and 
vice versa. Therefore, it could be argued that a detailed 
modelling of the leads is critical for the understanding of 
transport in graphitic systems. 

The objective of this paper is to show that to the 
contrary, transport through a sufficiently large and only 
weakly doped graphene sample often does not critically 
depend on most details of the leads. For the particular 
case of quantum-wire and graphitic leads, this insensitiv- 
ity manifests itself in an explicit duality, which can be 
formulated in terms of a shift of the gate potentials that 
control the charge density in the leads. This duality origi- 
nates in a mode selection mechanism, which dictates that 
the transport in extended, only weakly doped samples of 
graphene is dominated by the small part of the mode 
space in the vicinity of the conical points. Using this 
universal mechanism, one is naturally led to formulate 
an effective contact model, which is parameterized by a 
single complex number [i. This effective model can be ap- 
plied to a much larger class of leads which are amenable 
to the mode-selection mechanism. The main requirement 
is that the leads are sufficiently wide, and provide a dense 
set of propagating modes with transverse wave number 
close to the conical points (leads which do not possess 
any propagating modes at these transverse wave num- 
bers exhibit a very large contact resistance, and hence 
do not provide a good electronic coupling). 

In striking contrast to universal transport through 
semiconductor quantum dots^, the effective contact 
model provides a unified description of leads which sup- 
port a different number of propagating modes (but share 
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the same width). Numerical computations confirm that 
the effective contact model provides an efficient and 
reliable tool to study the transport through extended 
weakly-doped graphitic systems. 

This paper is organized as follows. Section HT1 formally 
poses the problem of the lead-sensitivity of transport in 
mesoscopic systems. In Section UTTT I calculate the con- 
ductance and Fano factor of a rectangular graphene sam- 
ple which is connected to quantum wires or graphitic 
leads. The rectangular graphene sample is held at the 
charge neutrality point, while the charge density in the 
leads is controlled by gate potentials. At a finite charge 
density in the leads, this deliver the duality of both con- 
tact models. In Section ITV| I describe the effective con- 
tact model which represents a large class of leads, en- 
compassing quantum wires and doped graphitic strips. 
Section [V] describes how the complex number \x could be 
extracted from the conductance and shot noise, and tests 
the applicability of the effective contact model via numer- 
ical computations for rectangular and circular graphene 
samples. Conclusions are presented in I VII The appen- 
dices contain technical details of the calculations. 
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II. FORMULATION OF THE PROBLEM 

At low bias, the phase-coherent transport properties of 
a mesoscopic system probed in a two-terminal geometry 
are encoded in the scattering matri x 13 i 14 



S 



r a 

t r' 



(1) 



which contains the transmission (reflection) amplitudes 
t, t' (r, r') of charge carriers incident from the source or 
drain contact, respectively. Two characteristic transport 
properties are the Landauer conductance 

G = (2e 2 /h)trth 



and the shot-noise Fano factor 



F = 1 



tritfttf 
tvtU 



(2) 



(3) 



In general, the transmission matrix t is determined by 
matching the propagating modes in the leads to the 
modes in the mesoscopic system, and hence depends on 
the properties of both of these constituents. This pa- 
per considers the case that the mesoscopic system is an 
extended, only weakly doped sample of graphene with a 
specified geometry, and asks the questions how the trans- 
port properties change when the sample is contacted by 
different types of leads (keeping the width and position 
of the leads fixed). 



III. DUALITY OF METALLIC AND 
GRAPHITIC CONTACTS 

It is instructive to first explore the potential relevance 
and eventual insensitivity to most details of the leads by 



FIG. 1: A graphene sample, modelled as a hexagonal lattice 
(lattice constant a), is attached to leads formed either by 
doped graphene (top) or by a quantum wire, modelled as a 
commensurably matched square lattice with lattice constant 
\/3a (bottom). The charge carrier density in the leads is 
controlled by a gate potential V g , while the central graphitic 
region is held close to the charge-neutrality point. The paper 
demonstrates that both leads are equivalent when the gate 
potential is suitably adjusted, and further extends this duality 
to a broad class of leads and contacts, which all share the 
properties of a simple effective contact model. 



calculating the conductance for a specific example, a rect- 
angular undoped graphene sample of width W = \/3aW 
and length C = 3aL, where W,L^>1 are integers (see 
Fig. [T]). The strip is connected to leads of the same 
width, either formed by doped graphene or by a quan- 
tum wire, which arc modelled in a tight-binding approach 
on a hexagonal or square lattice, respectively. The elec- 
tronic density in the leads is controlled by a gate poten- 
tial, inducing an on-site potential energy denoted by V g . 
For -f/W <§; \V g \ <C 7 the transport of this system con- 
nected to graphitic leads has been investigated earlier in 
the framework of the Dirac equation^ It was found that 
the conductivity a = = is of order of the exper- 
imentally observed value, while the Fano factor F = 1/3 
coincides with the universal value of a disordered quan- 
tum wire. 

A unified description of square- and hexagonal-lattice 
leads can be achieved when the hexagonal sublattices A, 
B, which differ in the orientation of the bonds, are in- 
dexed by two numbers I, m which are either both integer 
or both half-integer (see Fig. [T|). On the square lattice 
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the indices / and m are both integer. For both types of 
lattice, the transverse wave function of the modes in the 
lead are then given by 



square-lattice leads 



hexagonal-lattice leads 



w + 1 w + r 



(4) 



where n — 1, 2, 3, . . . , W is the mode index. Associated 
with each transverse mode are two extended Bloch waves 
with longitudinal wave numbers ±fe„, which are real for 
propagating modes and complex for evanescent modes. 
The wave numbers are fixed by the dispersion relation, 
which for the hexagonal lattice is given by 



Vg = vVfn,k„fn ,-*:„, r? = sgnV s = ±1, (5a) 



fn,k n = 7 + 2 7 e l3fc " a / 2 ■ 



2(W+1) 
while for the square lattice 

V g = 27 cos(\/3afc n ) + 27 cos 



W + l 



(5b) 



(G) 



Because of the different dispersion relations, the sets of 
transverse-mode indices n supporting propagating modes 
on the two types of lattice are not the same. Hence, for 
ordinary mesoscopic systems (such as a semiconducting 
quantum dot, a quantum wire, a metallic nanoparticle 
or a hybrid structure with superconducting or ferromag- 
netic properties) one would obtain distinctively different 
transport properties when quantum wire leads would be 
replaced by graphitic leads. 

In order to obtain the transmission amplitudes t n >n 
when these leads are coupled to the undoped rectangu- 
lar graphene sample, the modes in the leads have to be 
matched to modes in the graphene strip, which follow 
from Eq. §5§ by setting V g — 0. For the present quasi-one 
dimensional geometry, the mode index n is conserved,- 
and a straightforward calculation (see Appendix [AJ de- 
livers 



Jn'n f 



where the moduli of 



2 , ( n nir 

— In 2 cos — ^ 

3a V 2(W+l) 



(7) 



(8) 



denote the decay constants of the modes in undoped 
graphene. The lead dependence of the transmission coef- 
ficients is encoded in the self energies /-t n ,± of the in- and 
outgoing modes, given by 



..(h) _ 
Mn,± — 

u (s) - 
H"n,± ~ 



- 7 e 



fn,±7jk„ 



yj fn,k„fn,-k„ 

±iV3k n 



(hexagonal lattice), (9a) 
(square lattice). (9b) 



The effective contact model developed in Section IIVI 
arises from the observation that the detailed lead de- 
pendence embodied in these numbers becomes irrelevant 




FIG. 2: Gate-voltage dependence of the conductivity a — 
(C/W)G (in units of e 2 /h) and the shot-noise Fano factor F of 
an undoped graphene strip of width W = 152 \/3a and length 
C — 30 a. The open data points are obtained by numerical 
computations for square-lattice leads (left) and hexagonal- 
lattice leads (right). The solid curves are the prediction of 
the effective contact model ()18|) with \i given by Eq. (|10[) 
or Eq. respectively. The dashed curves are the analyti- 

cal expressions (|12I13I14I15[) obtained within the saddle-point 
approximation. 



when the leads are coupled to a sufficiently large graphitic 
region. According to Eq. ([5]), all modes in this region 
decay rapidly with the exception of modes with index 
71 2{W + l)/3 = h, which have transverse wave num- 
bers in the vicinity of the conical points. For a sufficiently 
large sample we hence only require to know the complex 
number /i = fifi,+ = A*s — ) which for a hexagonal-lattice 
lead is given by 



1 



V W (V g ) = --[V g+ix /4 1 i-V2 
while for a square-lattice lead 



(10) 



(V g ) = -\{v 9 + 7 + * V^ 2 - (V g + 7) 2 ) • (11) 



Before we describe the general consequences of this 
observation we first proceed to explore the consequences 
for the rectangular graphene sample. To make contact to 
Ref. d, let us further assume that the strip is very wide, 
W ^> C » a. The conductance can then be calculated 
in a saddle-point approximation, which is described in 
Appendix [B] For hexagonal-lattice leads this gives 



G = 



4e 2 W 



v 2 



irh C 



V a 



arcsm 



— , (12) 
2 7 ' V ' 



which for V g — > recovers the result G 
from the Dirac equation^ For square-lattice leads one 



^t-^t- derived 
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finds 

G = ^V^-^ + ^anrin^+I. (13) 
A similar calculation yields the Fano factor 

2 7 2 JA^ - V* 

F = — 1 (14) 

V* 2^ g arcsin(Vy 2 7 ) v ' 

for the hexagonal-lattice leads, and 

p = _jf V^ 2 -(v g +i) 2 , , 

(V g + -f) 2 2(l/ g +7)arcsin[(l/ ff +7)/2 7 ] 1 ' 

for the square-lattice leads. For a numerical validation 
of these results see Section [Vj 

Equations (|12|14|) for the hexagonal-lattice leads and 
(|13|15|) for the square-lattice leads coincide when the gate 
potential is shifted by 7. In particular, for the square- 
lattice leads the values G — of the conductance 
and F = 1/3 for the Fano factor are now recovered for 
V g — ► —7. This is a direct consequence of the relation 

^(V g )=^(V g+1 ) (16) 

between the characteristic self energies /i, Eqs. (fTT))) and 
(fTTj) . which describe the coupling of the leads to the con- 
ical points. This relation does not hold for the self en- 
ergies ([9]) away from the conical point, but this is not 
detected in the transport as long as the leads support 
sufficiently many propagating modes around n — n. Sur- 
prisingly, Eq. (fTB"]) entails that the point of half-filling of 
the hexagonal-lattice leads (the charge-neutrality point 
V g = 0) corresponds to the point of three-quarter filling 
in the square-lattice leads. 

IV. EFFECTIVE CONTACT MODEL 

In order to elucidate the generality of this duality be- 
tween hexagonal- and square-lattice leads, we now lib- 
erate ourselves from the strip geometry of the graphene 
sample, hence, consider samples of more arbitrary geom- 
etry in which different transverse modes can be mixed 
by the transport. In this case the matrix of transmission 
amplitudes t n '„ is no longer diagonal, and in general can 
be obtained from the Fisher-Lee relation^ 

S = -l + iVv¥v^(E-H-'E)- 1 r^Vv. (17) 

The leads are represented by their self energy £ = 
V^l diag [fx nt+ )^V\ where V is a coupling matrix of 
the leads to the contact region. The diagonal matrix 
v = diag (— 2 Im/i„ !+ ) contains factors proportional to 
the propagation velocity of the modes in the leads. The 
matrix '5 now accounts for the transverse modes in all 
leads. 



According to Eq. ([9]) , the modes in the hexagonal- and 
square-lattice leads in general have different propagation 
velocities and self-energies. As a consequence, when these 
leads are connected to an arbitrary system, described 
by the internal Hamiltonian H, the resulting transport 
properties in general will differ. When H represents a 
sufficiently large, weakly doped graphitic system, how- 
ever, the localization of the Fermi surface near the coni- 
cal points guarantees that only the lead modes coupling 
to this part of the Brillouin zone will contribute signifi- 
cantly to the transport. This also applies to other types 
of leads, different from the two specific cases considered 
so far, and then yields the universality among contact 
models advertised in the introduction. 

Assuming that the leads support a large number of 
propagating modes whose properties depend smoothly 
on the quasi-continuous mode index n, we hence can 
equip all modes in the lead with the same constant 
/i„ !+ — * /i = fin,+ characteristic of the modes coupling 
to the conical points. Using that the transverse mode 
profiles form an orthogonal set, = 1, the self energy 
then simplifies to £ = f/PV'. For optimally matched 
leads, V = P is just a projector of the internal system 
space onto the contact region. A unitary transformation 
S — > ^fS^ (which does not affect the conductance and 
the Fano factor) then results in a simplified Fisher-Lee 
relation for the scattering matrix, 

S = -1 - (2iImfi)P T (E - H - iiPP T )- 1 P. (18) 

This form of the scattering matrix defines an effective 
contact model, which is parameterized by a single com- 
plex number \i. 

For hexagonal-lattice or square-lattice leads, fi is given 
by Eq. (fTU]) or Eq. (fTTj) . respectively. However, our ar- 
gumentation entails that the simple effective model Eq. 
(fT5j) can be applied to a much larger class of leads, which 
share the same mode-selection principle when coupled to 
graphene. (An example would be a square lattice with 
a reduced lattice constant ^/3a/q, where q is an inte- 
ger.) Within this class of leads, [i is still given by the 
self energy Depending on the contacts, the ma- 

trix W' may not be uniform over the contact region, 
and also may account for a contact tunnel-barrier re- 
sistance. Assuming that the coupling strengths p n = 
{^\ l VV^ n )/{^iPP T ^ n ) (where is the column vec- 
tor associated to the nth transverse mode) do not depend 
strongly on the mode index, the effective contact model 
(fT5j) still holds with /i = pnfin y+ J£~ Ideal coupling into 
the graphitic modes near the conical points is realized 
for a self energy /1 = —27. 



V. APPLICATIONS 

This section describes some consequences of the effec- 
tive contact model (|18p . and compares its predictions to 
the results of numerical computations. 
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FIG. 3: Gate-voltage dependence of the conductance G (in 
units of e 2 /h) and the shot-noise Fano factor F for a weakly 
doped circular graphene sample (Ef = 7/10) of radius TZ = 
100 v3o, connected to leads of width W = 60 \/3ci- The data 
points are obtained by numerical computations for square- 
lattice leads (left) and hexagonal-lattice leads (right). The 
curves are the predictions of the effective contact model (|18p 
with /j, given by Eq. (flO|) or Eq. (fTTj) . respectively. 



Transport in a wide, long graphitic strip and 
the reconstruction of /1 



Within the effective contact model (fla|) . the conduc- 
tance of a wide and long graphene strip (as considered in 
Section [TTT]) is given by 



„ 4e 2 W Im a . Re n 

G = —r-r-^ arcsm— -, 

irh L tie /j, \fi\ 



while the Fano factor is given by 

iMl 2 Im, 



2(Re/i) 2 2Re/zarcsin(Re/u/|/x|) 



(19) 



(20) 



(for the derivation see Appendix[B]). These formulae gen- 
eralize Eqs. (|12|) - (fT5|) . which are recovered when /i is 
taken from Eq. (fTTj)) (for hexagonal-lattice leads) or Eq. 
(fTTj) (for square-lattice leads). For other types of leads, 
the two transport characteristics deliver two independent 
numbers which can be used to infer the complex number 
/i. 



B. Numerical results 

The validity of the effective contact model can be as- 
serted by numerical computations. The gate-voltage de- 
pendence of the conductivity a and the Fano factor F 
for a sample of width W = 152 a and length £ = 30 a 
is shown in Fig. f2J Results of numerical computations 
obtained by the method of recursive Green's functions 17 



are compared to the theoretical prediction of the effective 
contact model resulting from Eq. (fTU|) or Eq. (fTTj) . as well 
as to the analytical expressions (fT2")) - (fT5|) . Good agree- 
ment is found in the range \V g + 7I < 27 for the square- 
lattice leads, and 7/W < | V^, | < 27 for the hexagonal- 
lattice leads, corresponding to the condition that the 
modes in the vicinity of n are propagating. 

As an additional example, Fig. [3] shows the conduc- 
tance and the Fano factor of a weakly doped circular 
graphitic region (Ep = 0.1 7), calculated for hexagonal- 
and square-lattice leads and within the effective contact 
model (fTB|) , where fi is again given by Eqs. (fTTJ|) . (fTTj) . 
With the exception of the region \V g \ f» 7/W, where the 
hexagonal-lattice lead does not support many propagat- 
ing modes, the predictions of the effective contact model 
again agree nicely with the results for the two types of 
leads. This shows that the effective contact model re- 
mains applicable for graphitic samples where modes are 
mixed in the transport (so that the transmission matrix 
is no longer diagonal). 



VI. SUMMARY 

This paper critically assesses how sensitively the trans- 
port through extended samples of weakly doped graphene 
depends on the details of the leads and contacts connect- 
ing the sample to the electronic reservoirs. The starting 
point is the observation that there exists a duality be- 
tween doped graphitic leads and quantum wires, which 
result in the same transport properties if a gate voltage 
is suitably adjusted. This duality is not based on sim- 
ilarities of the two type of wires (which would result in 
very distinct transport properties if attached to a conven- 
tional material), but is rooted in a unique mode selection 
mechanism which originates from the conical points of 
undoped graphene. Maybe most strikingly, the duality 
holds even though both types of leads support a different 
number of propagating modes. 

Since the mode selection mechanism is a universal 
property of weakly doped graphene, the duality described 
above can be generalized into an effective contact model, 
which applies to a broad class of leads. The effective 
model is parameterized by a single complex number /z, 
which can be determined by a measurement of the con- 
ductance and the shot noise of a rectangular undoped 
graphene strip. 

In the present paper, this effective model has been de- 
rived using the technical requirement that the leads pro- 
vide a densely spaced set of propagating modes which 
couple smoothly to the conical points. This requirement 
is fulfilled for the typical leads considered in past and 
present mesoscopic-transport studies, and especially, for 
ballistic, wide contacts which provide a good electronic 
coupling to the graphitic sample (so that a transport 
measurement probes the sample, and not merely the con- 
tact resistance). 

One potential mechanism to violate the specific tech- 



6 



(a) 



(b) 





FIG. 4: (a) Brillouin zone for a square lattice with lattice 
constant y/3a. (b) Hexagonal Brillouin zone for a honeycomb 
lattice with lattice constant a (dashed) , and its rearrangement 
into a rectangular Brillouin zones (solid). The numbered tri- 
angles denote regions which are transported with reciprocal 
lattice vectors. The solid dots denote the K points. 



nical assumptions used to derive the effective contact 
model is strong interface disorder, so that the propagat- 
ing modes in the leads are mixed by the interface. For 
an ordinary mesoscopic system, such a diffusive interface 
changes the specific transport properties decisively (these 
changes cannot be described by a simple tunnel barrier). 
For weakly doped graphene, it is conceivable that even 
such a drastic modification of the contacts does not com- 
pletely violate the effective contact model. In a transfer 
description of the contact, the crossing of the diffusive 
interface from the graphitic sample to the lead results in 
a translation of the transverse mode profile of the modes 
close to the conical points into a random superposition 
of the modes in the lead. It is conceivable that self- 
averaging in these random superpositions merely renor- 
malizes the effective self-energy /j, into the average self- 
energy of the propagating modes in the lead. This ques- 
tion goes beyond the scope of the techniques used in the 
present paper, and is left for further consideration. 
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APPENDIX A: MODE MATCHING FOR AN 
UNDOPED RECTANGULAR GRAPHENE STRIP 

This appendix provides details of the derivation of the 
transmission eigenvalues (0 for an undoped rectangu- 
lar graphitic strip coupled to square-lattice or hexagonal- 
lattice leads. 

A unified description of both cases in real space is fa- 
cilitated by the discrete coordinate system (I, to), shown 



in Fig. [TJ Analogously, Fig. [4] shows a convenient choice 
of the Brillouin zones which facilitates a unified descrip- 
tion in k space. For the square lattice with lattice con- 
stant ^/3a, the Brillouin zone is given by the standard 
square |fc x |, \ k y \ < Using the periodicity of the band 
structure, the conventional hexagonal Brillouin zone for 
the honeycomb lattice is rearranged into the rectangle 
\k x \ < |f = Qx, \k y \ < ^ ee \q yi where \k x \ = q x 



\k y \ = q y denotes the position of the conical points. With 
this choice, both Brillouin zones have the same extent in 
k y direction. Moreover, the boundary conditions of the 
tight-binding model select the same transverse wavenum- 
bers k v n = ,„ 7m / - . The transverse mode profile ^ nm is 
then uniform throughout the system, and given by Eq. 
©■ 

For fixed transverse wavenumber k y ^ n , the longitudinal 
wavenumbers k n follow from the dispersion relations ([5|) 
and ([5]), respectively. For each transverse wavenumber, 
there are two solutions of opposite sign. 

In the undoped central graphene regions, the longi- 
tudinal wavenumbers are complex, and the modes are 
evanescent with decay constant Im k n — K n given in Eq. 
(JH). The wave functions <pA , 4>b on the A and B sites are 
determined by the Schrodinger equation, which delivers 



4% to) 



-3alK n ~\-27Til 



3alK, n -\-2nil 



(Ala) 
(Alb) 



respectively, where a and /3 are constants. 

The Landauer approach requires to match these wave- 
functions to the propagating modes in the lead. These 
modes are characterized by a real longitudinal wavenum- 
bers ±fc„, where for definiteness k n > is taken to be 
the solution in the right half of the Brillouin zone. 

On the hexagonal lattice, modes with a positive longi- 
tudinal wavenumber k n > are right-propagating when 
sgnV^ = rj > 0, (so that the Fermi energy lies below the 
conical point), while these modes are left-propagating for 
t] < (since the band structure above the conical point 
has the opposite curvature). For a particle incident from 
the left, the solution of the Schrodinger equation in the 
left lead takes the form 

+ Cr ml r ? (/„,_ t , fe „//„_ I)fc J 1 / 4 e - 34 " fe " ai * n „(A2a) 

+ Cr nn (fn,- n k n /UmJ 1/4 e- 34?,fc " a ^„ m , (A2b) 
while in the right lead it is given by 

Ct nn Tj(^f rLj — r jf in I fn,r]k n ) ^ C " ' ^ ^ ^ n7Tl 



(A2c) 



<j> { £ L (l,m) = Ct nn (fn^ k Jfn,- vk J 1/4 e 3 ^< l - L ^ nm . 

(A2d) 

Here C is an arbitrary coefficient, while r nn and t nn are 
the reflection and transmission coefficient, respectively. 
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On the square lattice, positive longitudinal wavenum- 
bers k n > are always associated with right-propagating 
modes. The wave function in the left lead can hence be 
written as 



($\l,m) = C e V3,/c„a(* + l/2% r 



+ Cr nn e-^ lk ^ l+1 ' 2 ^ nm , (A3a) 



while in the right lead 



<jy^{l,m) = Ct nn e 



y3ik n a(l— L- 



1/2)^ 



(A3b) 



Using the fact that the wavefunctions I|A1|I . (|A2|) . jA3j, 
all satisfy the Schrodinger equation when the graphitic 
strip or the leads would be formally extended beyond the 
interface, the matching conditions take the simple form 
of continuity requirements 



4 n) (0,m) = 
4 n) (0,m) = 
4 n) (L,m) 



ft(0, m), 



4>b(L, m) = 4>^ R (L,m) 
for the hexagonal-lattice leads, and 



(0,m), 



4 n) (0,m) = ^"V L.,„). 
4 n) (0,m) 

4 n) ' 



\L,m) 
\L,m) 



4 n) (i,m), 

4 n) ' 



^(I + l.m) 



(A4a) 
(A4b) 
(A4c) 
(A4d) 

(A5a) 
(A5b) 
(A5c) 
(A5d) 



for the square-lattice leads. 

The resulting linear systems of equations deliver the 
transmission coefficient t nn in the form ([7]), where for 
each lattice the self energy is defined in Eq. |9]). 



leads follow when these expressions are evaluated with 
(fTUf and (fTTj) . respectively. 

Let us start with a detailed account for the conduc- 
tance, which in general is related to the transmission 
coefficients t nm via the Landauer formula f2|). The 
transmission matrix (J7|) of the graphene strip is diag- 
onal, hence G — (2e 2 /h) J2 n \t nn \ 2 . Since the strip is 
wide, W ^> C a, the transverse wavenumbers are 
closely spaced and the transmission coefficients ([7} de- 
pend quasi-continuously on the mode index n, so that 
the sum can be replaced by an integral, 



G 



2e2 
h 



f^n,-\- Mr; 



dn. 



(Bl) 



For C 3> a, all transmission coefficients decay rapidly 
with the exception of modes with index n » 2(W^+l)/3 = 
n, which have transverse wave numbers in the vicinity of 
the conical points and hence a small decay coefficient 
k„, given in Eq. ([8]). Hence the integrand in Eq. (|B 1 1) 
has a pronounced maximum at n m fi. Away from this 
maximum, for n = n + (5, the integrand decays rapidly 
because of the finite decay constant K n ~ #7r/W, and 
attains very small values on a scale on which the self- 
energies Hn,+ = H>n - barely change. One hence can ap- 
proximate the self energy by the value /J, n .+ ~ = A* 
at the conical point. Furthermore, the integration limits 
can be extended to ±oo. The resulting integral 



G- 2e2 



2(Im/i) 2 



dS 



, \n\ 2 cosh(2(57r£/>V) - (Re^) 2 + (Im/i) 2 

(B2) 

can be calculated exactly, and results in Eq. (fl9|) . 

For the Fano factor ([3]) one is led to calculate a similar 
integral I = J |inn| 4 dn, which can be evaluated using ex- 
actly the same steps as for the conductance. The integral 
is then approximated as 



APPENDIX B: CONDUCTANCE AND FANO 
FACTOR FOR WIDE AND LONG GRAPHENE 
STRIPS 

This appendix describes the saddle-point approxima- 
tion which is used to obtain the conductance (fTl?)) and the 
Fano factor (|20[) of an undoped graphene strip of width 
W > C > a. The conductance JEJ), ([13]) and the Fano 
factor (fT4|) . (fl"5|) for hexagonal-lattice and square-lattice 



/ = 



4(Im/x) 4 



[M 2 cosh(2<$7r£/W) - (Rc^) 2 + (Im^) 2 ] 2 



which again can be evaluated exactly, 



dS, 
(B3) 



W (Im/i) 2 W Im/4(Im/z) 2 - (Re/i) 2 ] . Refj, 
I = — - — — — — arcsm -; — — . 



ir£ (Re fi) 2 nC 



(Re fif 



m 

(B4) 



The final result for the Fano factor is given in Eq. (|20|) 
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